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I present a brief summary of the most important updates in the MRST determina- 
tions of the NLO and NNLO parton distributions. I concentrate on uncertainties 
due to input assumptions and possible theoretical corrections, particularly isospin 
violation, QED corrections, NNLO and corrections at low Q 2 and small x. 

1 Introduction 

A knowledge of the partonic structure of the proton is an essential ingredient in the 
analysis of hard scattering data from pp or pp or ep high energy collisions. Much at- 
tention has recently been devoted to obtaining reliable uncertainties on the parton 
distributions obtained from global or semi-global fits. One obvious uncertainty is 
due to the systematic and statistical errors of the data, i.e the experimental errors. 
They have been estimated by several groups PP-^O], working within a NLO (and 
sometimes NNLO) framework using a variety of different procedures. The general 
conclusion is that in all approaches the uncertainties are of the order of 2 — 5%, 
except for certain partons and regions of phase space, where the uncertainty can 
be bigger, particularly the high-x gluons and very high-x down quarks. However, 
when the partons and predictions for physical quantities using different parton sets 
are compared, the deviations can be significantly greater than the quoted errors. 
A recent summary can be found in jll) . This unfortunate discrepancy is likely to 
be due to different assumptions made by different groups when performing their 
fitting procedures, such as different data sets used, parameterizations for input sets, 
choices of heavy flavour prescriptions, ordering of perturbative series, etc. However, 
the different results for partons obtained when making different assumptions also 
suggest that the standard fixed order in the perturbative expansion is not com- 
pletely adequate to describe the data, with different individual data sets making 
different and conflicting demands on the partons and on as- Indeed, MRST do see 
certain difficulties in obtaining the best possible fit to the data with the current 
theoretical treatment - the HERA and NMC data in the range x ~ 0.01 — 0.001 
are too steep in Q 2 ; the gluon is negative (or at best valencelike) at small x and 
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Q 2 ; and the Tevatron jet data prefer a different shape gluon from the DIS data. 
Hence, in order to understand the uncertainties in parton distributions, it is vital to 
consider theoretical corrections. These include the possibility of isospin violation; 
s(x) ^ s(x); higher orders (NNLO); QED (comparable to NNLO ? i.e. a 3 s ~ a); 
large x «ln 2n_1 (l - a;)); low Q 2 (higher twist); small x « ln n_1 (l/x)). The 
uncertainties associated with many of these possible sources have recently been 
considered ^21- Here I will highlight some of the main conclusions and some addi- 
tional results which postdate this analysis. 

One of the issues we have examined is the possibility of isospin violation. This 
is particularly relevant because of the "NuTeV anomaly" NuTeV measure, 

with a 3 — ex discrepancy, 

where [SU V ] = [U$] - [£>£] and [SD V ] = [D$] - [U?], [U$] = J^x U P{x)dx etc.. We 
parameterize the isospin violation by 

u*(x) = d»(x) + K/(a;), dP(x) = <(*) - Kf(x), 

where f(x) is a number conserving function with valence like behaviour at large and 
small x 52 ■ The variation of the x 2 f° r the global fit is shown in Fig. 1 , and the 
best fit leads to k = — 0.2 (which agrees very well with a prediction in Jl]), and to a 
reduction of the NuTeV anomaly to ~ 1.5 — it, i.e. Asin 2 &w ~ —0.002. However, a 
wide variation in k is possible, and the anomaly may be reduced further or enlarged. 
Other issues to which the anomaly is sensitive are the strangeness asymmetry in the 
proton recently investigated in ^S] and electroweak corrections I n each case 
the uncertainties may be large as the anomaly, though there is some disagreement 
regarding the strangeness asymmetry. 
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Figure 1. Variation of global x 2 with the value of k 

We have also recently examined the effects due to QED. The corrections to 
DIS include those in Fig. 2 and lead to mass singularities when 7 || q. The QED 
collincar singularities are universal and can be absorbed into pdfs, exactly as for 
QCD collinear singularities, leaving finite (as m q — > 0) O(a) QED corrections in 
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coefficient functions, and QED-improved DGLAP equations. The effect on the 
gluon distribution is negligible as is that on the quark distributions at small x 
where the gluon contribution dominates DGLAP evolution. Only at large x is the 
effect noticeable (0(1%)) as first noted in [TJI, an d more recently verified in |18|. 
and it is equivalent to a slight shift in as- Aas(Mj) ~ +0.0003. This is much 
smaller than many sources of uncertainty. 





Figure 2. Diagrams for QED corrections to DIS 

However, these QED effects, although small, lead to isospin violation. Uy(x) 
quarks radiate more photons than dy(x) quarks. To a rough approximation, the 
photon distribution should be 



7(x, 



z!7r J x v v 



So there is more photon momentum in the proton than in the neutron due to high-x 
up quarks radiating more than high- a; down quarks. Momentum conservation leads 
to Uy{x) < d v (x) at high x. Hence, [SU V ] < as required by the NuTeV anomaly 
and as found by our best fit above. Our estimates imply that the QED effects lead 
to ~ 3/4 of the isospin violation observed by our best fit. Details of this and parton 
distributions involving QED corrections will soon appear |19|. 

Much more generally important are NNLO corrections. The coefficient func- 
tions for F2(x,Q 2 ) are known at NNLO [201 . Singular limits x — > and large Np 
limits are known for the NNLO splitting functions |23 _ E3 as wei l as limited mo- 
ments and approximate NNLO splitting functions were devised and used 
by MRST g£| (and Alekhin 0) in NNLO fits. Compared to NLO they improved 
the quality of the fit very slightly, and reduced as to ~ 0.1155. The predictions 
lead to reasonable stability order by order for (quark-dominated) W and Z cross- 
sections. However, this fairly good convergence is largely guaranteed because the 
quarks are fit directly to data, and it is much worse for gluon-dominated quanti- 
ties e.g. Fl(x,Q 2 ) or Higgs production. Very recently, the full calculation of the 
NNLO splitting functions was completed We have confirmed that these alter 
the parton distributions obtained from the previous best estimates by only a very 
small amount. Full NNLO partons will appear soon, but also require the complete 
incorporation of NNLO coefficient functions for as many processes as possible, and 
a rigorous treatment of heavy flavour thresholds at NNLO [25J . For predictions it 
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will be useful to have the exact NNLO (0(ctg)) longitudinal coefficient functions, 
and these, along with the 0(a? s ) (NNNLO) coefficient functions for Fi(x, Q 2 ) will 
soon be complete. We have estimated the most likely value, along with uncertain- 
ties, from the limited information, as had previously been done for the splitting 
functions. The central value and the two "extremes" for the dominant contribu- 
tion, i.e. the part of C\ g which is cx Nf are respectively (the common factor of 
^/(oW^tt)) 3 is taken out and L = ln(l/x), L x = ln(l - x)) 

L 1212 9 , 

410 1685 - 54161s - 47131Li - 11496L? - 867L?, 

x x 

410 + 20298 - 106974x - 67527Li - 16336L 2 - 1241L 3 - 14116X, 

x x 
T on 

410 — + 30895a; 3 + 6220Li - 102L? - 66.6L? + 1378L 2 . 

x x 

The failure of the NNLO corrections to solve completely the problems with the 
quality of the fit leads us to try an alternative approach. In order to investigate 
the real quality of the fit and the regions with problems we vary kinematic cuts on 
data, i.e. change W 2 ut , Q 2 ut and x cutl re-fit and see if the quality of the fit to the 
remaining data improves and/or the input parameters change dramatically. This 
is continued until the quality of the fit and the partons stabilize. This is discussed 
in detail in [T2]. Briefly, our W 2 cut of 12.5GeV 2 seems reasonable, whereas the fit 
improves gradually until Q 2 ut ~ lOGeV 2 and x cu t ~ 0.005.° This is true at both 
NLO and NNLO, but the modification of the partons with changing cuts is much 
greater at NLO. For both NLO and NNLO we have defined "conservative" partons 
12 which arc obtained using only the data within the conservative cuts. 

In order to investigate the possible theoretical improvements which could modify 
the partons at small x and low Q 2 (particularly at high x) we have tried including 
known ln(l— x) resummations, saturation corrections, and purely phenomenological 
higher twist and ln(l/a;) corrections. At high x, higher orders in the coefficient 
functions and higher twist contributions can improve the fit, but are difficult to 
distinguish. At low x, there is no evidence for phenomenological higher twist. 
Saturation corrections [SDJ do not seem to help the NLO b or NNLO fits (they 
would help only at LO). There is a distinct improvement in the fit quality when 
(momentum conserving) ln(l/a;) corrections are added to the splitting functions, 
increasing the small- a; evolution. Hence, there is some reason for uncertainty about 
the correct approach, and hence the correct predictions at very small x. Indeed, a 
comparison of the LO, NLO, approximate NNLO and a resummed prediction [22] 
for Fl(x,Q 2 ) is shown in Fig. 3, exhibiting the level of uncertainty. Also shown 
are potential data which could be obtained by low energy running at HERA. It is 
clear that such data would play a huge role in enabling us to understand and be 
confident of predictions involving small x partons. 

In summary, as well as the uncertainties on partons due to experimental errors, 
there are uncertainties on partons, and related physical quantities, due to assump- 
tions in the fit and having to work within an incomplete theoretical framework. 



a We note that CTEQ do not seem to reproduce these conclusions 1291 . 

b A recent joint analysis with diffractive data suggests some slight improvement 1311 . 
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Some of the uncertainties are small, but in some cases, e.g. QED corrections or 
s(x) 7^ s(x), they can be important for specific quantities. Others, associated with 
low Q 2 and small x, may be potentially large, and more theoretical work is needed, 
along with specifically targeted data to guide and verify this work. 



F L LO , NLO, NNLO and resummed - Simulation of Low Ep HI Data 



Q =5 GeV" 



NLO fit 
NNLO fit 
resum fit 
',LO fit 



a 

of, 




10 ' 10 10 ~ 10 10 1 



Q-=20 GeV" 




I; 1 " 


Q 2 =10 GeV 2 


_ ""if 




1 mi 


J 1 1 







5 io 4 io : 


10 2 10 1 


{"in " i 
i \ \ '• q2 


1 1 

=40 GeV 2 


; % 

i i 





10 10 10 ~ 10 10 1 10 10 10 ~ 10 10 



Figure 3. Predictions for Fl(x, Q 2 ) at LO, NLO, NNLO and with a resummed calculation 
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